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ABSTRACT 

Extinction  distance  measurements  provide  a  method  of  testing 
the  predictions  of  the  dynamical  theory  of  electron  diffraction# 

Using  silicon  as  the  specimen  material,  we  have  made  measure¬ 
ments  of  the  variation  of  the  (111)  and  (220)  extinction  distances  with 
the  angular  deviation  of  these  reflections  from  their  Bragg  orientations# 
The  (220)  measurements  showed  good  agreement  with  both  the  two-beam 
and,  taking  into  account  systematic  reflections,  multi-beam  approxima¬ 
tions  of  the  dynamical  theory#  However,  the  (111)  results  were  in  poor 
agreement  with  the  two-beam  approximation  and  only  in  qualitative  agree¬ 
ment  with  the  multi-beam  approximation*  Better  agreement  was  found  with 
a  multi-beam  approximation  in  which  an  experimental  value  of  2 #83  A 
for  the  (111)  electron  scattering  factor  was  used.  This  value  was 
calculated  from  the  value  of  840  £  40A  measured  here  for  the  (111) 
extinction  distance  at  152  kV.  These  latter  calculations  are  in  good 
agreement  with  recent  measurements  of  the  corresponding  X-ray  scatter¬ 
ing  factor# 
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CHAPTER  I 
INTRODUCTION 


The  interpretation  of  the  results  of  electron  diffraction 

experiments  depends  on  the  type  of  interaction  the  incident  electrons 

are  assumed  to  have  with  the  target  material®  Until  the  perfection 

of  the  electron  microscope"^  electron  diffraction  results  were  usually 

interpreted  in  terms  of  the  kinematical  theory  of  electron  diffrac® 

tion  which  is  similar  in  its  approach  to  the  theory  commonly  used  in 

X~ray  diffraction®  However,  electron  microscope  images  of  wedge® 

2«? 

shaped  crystalline  materials  obtained  in  the  early  1940‘s  had  in® 
tensity  variations  which  could  only  partially  be  explained  by  this 
theory.  These  variations,  known  as  thickness  extinction  contours 
were  similar  to  those  shown  in  Figure  1.  The  kinematical  theory  was 
able  to  account  for  these  observations  by  predicting  that  the  dif® 
fracted  beam  would  vary  periodically  with  depth  in  the  crystal.  As 
shown  in  Figure  2,  this  periodicity  would  result  in  the  presence  of 
extinction  contours  in  dark  field  electron  micrographs  of  wedge-¬ 
shaped.  crystals.  However,  th©  kinematical  theory  also  predicts  that 
the  period  of  this  variation  or  the  extinction  distance  becomes  in® 
finite  as  the  diffracted  beam  approaches  the  Bragg  condition®  This 
was  in  disagreement  with  experimental  observations  "  that  the  ex¬ 
tinction  distance  remained,  finite®  Because  of  the  failure  of  the 
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Figure  1 

Dark  field  image  of  a  silicon  crystal  showing  the  presence  of  extinc¬ 
tion  contours  in  wedge-shaped  regions*  x  20,000 
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Figure  2 

A  schematic  illustration  showing  the  origin  of  extinction  contours 

(a)  Cross  section  of  a  wedge-shaped  crystal  showing  the  depth  perio¬ 
dicity  of  Igt  the  intensity  of  a  diffracted  beam,  6^  is  the  extinc¬ 
tion  distance#  The  arrow  points  in  the  direction  of  the  incident  beam# 

(b)  The  dark  field  image  of  the  crystal  with  extinction  contours  at 
positions  E 


FIGURE  2 


kinematical  theory  near  the  Bragg  orientation,  experimenters  turned  to 
the  dynamical  theory  of  electron  diffraction  in  order  to  interpret 
their  results.. 

Q 

The  dynamical  theory  had  been  proposed  by  Bethe  in  1929  but 
had  been  largely  ignored  in  favour  of  the  simple?  kinematical  theory* 

In  his  theory  Bethe  gave  a  wave  mechanical  treatment  of  the  motion  of 
a  high  energy  electron  in  the  periodic  field  of  a  crystal*  This  treat- 
ment  resulted  in  an  infinite  set  of  equations  relating  the  amplitude 
coefficients  of  the  different  diffracted  beams  and,  thus,  no  general 
expression  for  the  intensity  of  an  image  could  be  obtained.  However, 
if  the  same  treatment  is  applied  to  a  finite  number  of  beams,  the  equa¬ 
tions  can  be  solved  numerically  to  give  the  intensities  of  the  beams. 

If  only  two  beams,  the  directly  transmitted  beam  and  one  strongly 
diffracted  beam,  are  considered,  the  equations  can  be  solved  analyt¬ 
ically,  Both  approximations  predict  that  the  intensity  of  a  strongly 
diffracted  beam  will  be  periodic  with  depth  in  the  crystal  even  at 

the  exact  Bragg  condition.  Thus,  the  dynamical  theory  can  explain, 

2-7 

at  least  qualitatively,  the  aforementioned  experimental  results. 

The  dynamical  theory  predicts,  further,  that  the  extinction 
distance  of  a  diffracted  beam  is  dependent  on  both  the  energy  of  the 
incident  electrons  and  A9,  the  angular  deviation  of  the  crystal  from 
the  exact  Bragg  orientation  for  that  reflection,  Relatively  little 
work,  however,  has  been  done  to  test  these  predictions  quantitatively. 
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q 

The  energy  dependence  has  been  examined  by  only  two  groups.  Hashimoto 
measured  the  variation  of  extinction  distance  in  aluminum  foils  with 
electron  energy  for  accelerating  voltages  from  50  to  300  kV  and  found 
that  the  experimental  results  agreed  very  well  with  the  theoretical 
two-beam  predictions.  However,  Dupouy^  et  al.  in  taking  similar  meas¬ 
urements  on  MgO  in  the  range  from  100  to  1000  kV  found  significant 

deviations  from  the  two-beam  predictions,  especially  at  the  higher 

11 

accelerating  voltages.  Goringe  et  al.  showed  that  these  deviations 
could  be  explained  by  taking  into  account  systematic  reflections  by 
means  of  the  multi-beam  dynamical  theory.  The  dependence  of  extinc¬ 
tion  distance  on  AO  has  been  investigated  by  only  one  group.  Uyeda 
12 

and  Nonoyama  measured  the  variation  of  extinction  distance  with 
angular  deviation  from  the  exact  Bragg  orientation  for  the  (200)  re¬ 
flection  in  MgO  wedge-shaped  crystals  and  found  good  agreement  with 
the  theory.  Their  results  were  obtained  at  random  A0*s  as  determined 
by  measurements  from  electron  diffraction  patterns.  However,  such 
determinations  are  of  only  moderate  accuracy,  a  fact  which  becomes 
very  important  when  small  A0*s  are  considered. 

Because  of  the  limited  amount  of  work  that  has  been  done, 
further  experiments  to  test  the  predictions  of  the  dynamical  theory 

J;  u 

were  proposed.  The  purpose  of  the  experiments  was  to^investigate  in 
a  systematic  and  accurate  manner  the  AO  dependence  of  extinction  dis¬ 
tance  and  to  compare  the  results  with  the  predictions  of  the  two-beam 


and  multi-beam  dynamical  theories  of  electron  diffraction*  It  is  a 


report  of  this  work  that  forms  the  basis  of  this  thesis. 
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CHAPTER  II 


DYNAMICAL  THEORY  OF  ELECTRON  DIFFRACTION 


The  dynamical  theory  of  electron  diffraction  as  first  pro™ 


8  13 

posed  by  Bethe  and  further  developed  by  MacGillavry  ^  and  Heiden- 

14 

reich  starts  with  the  Schroedinger  equation  for  an  electron  in  a 
crystal  potential  V(r)<*  This  equation  can  be  written 


V2>Kr)  + 


E  +  V(r)  Wr)  =  0, 


where  V(r)  is  the  wave  function  of  the  electron*  r  is  the  position 
vector,  m  and  e  are  the  rest  mass  and  charge  of  the  electron,  E  is 
the  potential  through  which  the  electron  was  accelerated  before  enter** 
ing  the  crystal  and  h  is  Planck's  constant*  Due  to  the  periodic  na¬ 
ture  of  the  crystal  potential,  V(r)  can  be  expressed  by  a  Fourier 
series 


where  v0  is  the  mean  inner  potential,  g  is  a  reciprocal  lattice  vec- 


means  that  the  term  g  =  0  is  excluded  from  the  summation.  Through  the 
use  of  Poisson's  equation  and  X-ray  scattering  factors,  the  Fourier 
coefficients  of  the  lattice  potential  can  be  found  in  terms  of  electron 
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scattering  factors. 


The  expression  is 


15 


4?.8? 


£ 

3 


fj(s) 


e-2rri  (g»Fj) 


volts. 


where  47.87  is  a  numerical  value  resulting  from  the  evaluation  of  phys- 
ical  constants,  £1  is  the  volume  of  the  unit  cell  in  A J  and  f  .(s)  is 

J 

the  electron  scattering  factor  in  A  of  the  j'th  atom  in  the  unit  cell. 
The  parameter  s  is  equal  to  4tt  sin  0/X  where  20  is  the  scattering  angle 
and  X  is  the  wavelength  of  the  incident  electrons. 

In  the  case  of  electrons  with  energies  of  the  order  of  those 
usually  used  in  electron  microscopes,  E  is  much  greater  than  V(r)  and 
th©  Schroedinger  equation  can  be  solved  using  a  perturbation  approach. 
The  solutions  of  the  equation  outside  the  crystal  where  V(r)  =  0  are 
plane  waves  of  the  form 

t|>(r)  ss  exp(2uiX«f) , 

where  the  magnitude  of  the  wave  vector  X  is  such  that 


2m 

The  effect  of  the  lattice  potential  is  to  give  solutions  of  the 
Schroedinger  equation  which  are  combinations  of  free  electron  plane 
waves o  Since  the  potential  is  periodic,  the  solutions  must  be  Bloch 
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functions,,  Thus,  they  are  of  the  form 

41(f)  =  |  C  (k)e2ni  (Vr) 

where  kg  s  k  +  g»  Ic  is  the  wave  vector  in  the  crystal  in  the  direction 
of  the  directly  transmitted  wave*  These  expressions  for  V(r)  and  vp(r) 
are  then  substituted  into  the  Schroedinger  equation  to  obtain 

e2ni(kg.f)  =  0( 

where  the  prime  on  the  second  summation  means  the  term  h  =  0  is  omit¬ 
ted,  Since  the  ^^(kg^r)  are  linearly  independent,  the  coefficients 
of  these  exponentials  must  all  equal  zero.  This  condition  results  in 
a  set  of  equations  of  the  form 

(K2  -  kg2)Cg(k)  +  f\cg_h(k)  =  0  C1) 

where  K  =  ~~  (E  +  v  )  and  Ue  =  ~  v_. 

h2  0  S  h2  g 

8 

This  set  of  equations,  called  the  "dispersion  equation"  by  Bethe  , 
gives  the  general  relations  between  the  amplitude  coefficients  C  ,  the 

o 

Fourier  coefficient  terms  U  ,  and  the  wave  vectors  of  the  electrons  k  • 

g  g 

It  is  from  these  relations  that  the  approximate  solutions  of  the 
dynamical  theory  are  developed. 


£ 

g 


g 


•cg(k) 


+  t[E 


v0]cg(k) 


2me 

h2 


h  Cg~h 


(k)v. 


« 
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A»  Two-Beam  Approximation 

In  the  two-beam  approximation  only  the  transmitted  beam  k 
and  one  diffracted  beam  kg  are  considered,,  In  this  case  the  set  of 
equations  reduces  to 

(K2  -  k?)C0(k)  +  U.gCg(k)  =  0 


(K2  -  kg2)Cg(k)  +  UgC0(k)  =  0 

Since  the  potential  energy  of  the  lattice  must  be  real 

ug  =  u-g* 

Also,  since  the  type  of  crystal  in  which  we  are  interested  is  centro- 
symmetric 

ug  =  u-g 

Therefore,  the  constant  Ug  is  real«  Now,  non-sero  solutions  exist 
for  CQ(k)  and  Cg(k)  only  if  the  determinant  formed  by  the  coefficients 
vanishes,  i.e*  if 

(K2  -  k2)(K2  -  kg2)  -  Ug2  =  0 

Since  K  and  k  are  nearly  equal  and  much  larger  than  g  in  magnitude, 
this  equation  can  be  rewritten  as 

U  ^ 

(k  -  K)(k  -  K)  = 

g  4K2 


< 


>  r  \ 


11 


The  values  of  k  satisfying  this  equation  lie  on  a  surface  in  k-space 

called  the  dispersion  surface  (see  Figure  3),  Physically,  the  theory 

predicts  that  a  degenerate  energy  level  corresponding  to  energy  iLiL, 

2m 

is  split  by  the  perturbing  potential  into  two  new  levels,  one  slightly 
below  and  one  slightly  above  the  old  level.  An  electron  in  the  crys¬ 
tal  with  wave  vector  K  is  thus  scattered  into  two  new  states  with 
wave  vectors  k^  and  k.^,  These  states  have  the  same  total  energy 
but  different  kinetic  energies. 

In  order  to  find  th©  amplitude  of  the  diffracted  beam  at  the 
exit  surface  of  the  crystal,  the  boundary  conditions  at  the  entrance 
and.  exit  surfaces  must  be  considered.  Since  the  energy  of  the  elec¬ 
trons  is  much  greater  than  th©  potential  energy  of  the  lattice,  re¬ 
flected  waves  from  the  two  surfaces  can  be  ignored..  The  boundary  con¬ 
ditions  then  reduce  to  simple  continuity  of  transmitted  and.  diffracted 
wave  amplitudes  at  th©  surfaces.  Thus,  the  amplitude  of  any  beam  just 
outside  a  surface  is  equal  to  the  sum,  taking  into  account  phase  dif¬ 
ferences,  of  the  amplitudes  of  the  Bloch  wave  components  in  the  direc¬ 
tion  of  that  beam  in  the  crystal* 

These  boundary  conditions  result  in  an  expression  for  the 
amplitude  of  the  diffracted  beam  which  is  proportional  to  the  differ¬ 
ence  of  two  exponentials.  If  the  surfaces  of  the  crystal  lie  approx¬ 
imately  parallel  to  th©  xy  plane  and  the  beam  traverses  it  in  the 
negative  z  direction,  this  difference  term  can  be  written 


. 
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Figure  3 

/T\  /p) 

The  dispersion  surface  in  the  two=beam  approximation*  '  and  D'  ' 
are  the  lower  and  upper  branches  of  the  surface  respectively*  In  this 
drawing  =  EH,  =  FE,  and  Ug/2K  =  AB  =  BC. 
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Sphere  of  radius  K 
centred  at  G 


Sphere  of  radius  K 
centred  at  O 


Reflecting  Sphere  of 
radius  K  centred  at  E 


FIGURE  3 
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m, 
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exp(2iT±  kg^z)  -  exp(2rri  kg^z) 
From  Figure  3  it  can  be  seen  that 


and 


kg^,  es  K  -  Y^^cos@g  ~  s  COS0g 

kg^^  a  K  +  Y^^COSQ  -  S  COS0 

0  o  o 


where  ~s  a  ~gt&*  The  general  sign  convention  for  s  is  that  s  is  nega<° 
tive  if  the  reciprocal  lattice  point  lies  outside  the  Ewald  sphere 
and  positive  if  it  lies  inside 0  In  the  case  of  electron  diffraction 


with  high  energy  electrons  K 


» 


g 


so  that  0  is  very  small  and 

O 


cosGg  can  be  set  equal  to  one*  From  geometrical  considerations  in 
Figure  3»  ©ne  can  show  that 


and 


XD„ 


s 

2 


+ 


s 

2 


Substituting  these  into  the  expressions  for  kg 


and  k we  obtain 

o 
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and 


(2)  _ 


=  K  + 


s 

«K- 

2 


When  these  expressions  for  k  ^  and  k  c '  are  substituted  into  the 

o  o 

exponential  difference  factor,  the  amplitude  is  found  to  be  propor= 
tional  to 


exp 


2tt i  (K 


s\_ 


sin 


nz( 


K2 


°r  Ig. 


the  intensity  of  the  diffracted  beam  is  proportional  to 


sin 


,u  g2 


At  the  exact  Bragg  condition  s  =  0  and  the  extinction  distance,  £,c, 
equals  K/U  *  C,  the  extinction  distance  when  s  is  not  equal  to  aero, 

o 

can  be  written  in  terms  of  in  the  form 


e. 


~  C0(l  + 


52K2' 


u, 


2 


or 


=  i/(i  + 


s2K2 

°g2 


* 

) 


This  last  expression  gives  the  variation  of  extinction  distance  with 


- 


' 


■ 
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angular  deviation  from  the  exact  Bragg  condition  as  predicted  by  the 
two-beam  dynamical  theory  of  electron  diffraction , 


B.  Multi-Beam  Approximation 

In  the  case  where  n  beams  are  taken  into  consideration  sev¬ 
eral  authors^”^  have  proposed  scattering  matrix  formulations  of  the 

17 

theory.  The  formulation  proposed  by  Howie  and.  Whelan  '  gives  the  most 
exact  numerical  solutions.  It  starts  at  (1),  the  dispersion  equation. 
If  the  determinant  of  the  coefficients  of  the  amplitudes  is  set  equal 
to  zero,  an  equation  of  degree  n  in  k  results.  The  solutions  to  this 
equation  define  n  dispersion  surfaces  in  k-space.  Then,  by  analogy 
with  the  two-beam  case,  for  the  i‘th  dispersion  surface 

K2  -  k^  =  -  2Ky^ 
and 

2 

K2  -  k^1)  =  2K(Sg  - 


Substitution  of  these  expressions  into  (1)  results  in  an  eigenvalue 
equation  of  the  form 


(2) 


where  is  a  column  vector  whose  elements,  Cg^,  are  the  wave  ampli¬ 
tudes  of  a  Block  wave,  is  th0  corresponding  value  of  y  and  Ais  a 
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matrix  whose  elements  are  A00  =  0,  Agg  s  sg  and  Agh  s  Ugc=h/2K  where 
g  ^  h. 

Consider  now  the  situation  shown  in  Figure  4.  Waves  in  the 
directions  of  the  n  beams  being  considered,  are  incident  upon  a  slab 
of  crystal  of  thickness  These  waves  can  be  represented  by  a  col¬ 

umn  vector  u  whose  components  are  the  amplitudes  of  the  various  waves • 
The  boundary  conditions  at  the  upper  surface  of  the  slab  result  in  a 
set  of  equations  of  the  form 


or,  in  matrix  notation 

ct  =  u, 

where  C  is  a  matrix  whose  columns  are  the  eigenvectors  of  (2)  and  Y  is 
a  column  vector  whose  elements  give  the  excitation  of  the  Bloch  waves 
on  the  various  branches  of  the  dispersion  surface.  ^  can  be  found 
from  the  equation 

T  ”1 

^  =  C  u  =  Cu  , 

since  C  is  an  orthogonal  matrix.  The  boundary  conditions  at  the  bottom 
surface  of  the  slab  are,  taking  into  account  phase  differences  of  waves 
from  different  dispersion  surfaces 

ugc  =  ?»{A^Cg^e:xp(2TTi  y^6z) 


I 


■ 
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Figure  4 


Waves  propagating  through  a  thin  slab  of  crystal 

(From  Howie  and  Whelan‘S) 


FIGURE  4 
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or,  in  matrix  form 

ua  as  Cy^', 

where  y  is  a  diagonal  matrix  whose  components  are  exp(2TTi  y^6z) „ 
Substituting  for  with  the  expression  we  obtained  above ,  we  have 

u»  =  Pu, 

,  -1 

where  JP  =  CyC  „  _P  is5  thus,  a  scattering  matrix  relating  the  ampli- 
tude  vectors  ua  at  the  bottom  surface  to  the  incident  wave  u  at  the  top 
of  the  crystal*  For  any  finite  number  of  beams,  the  components  of  the 
scattering  matrix  P  can  be  computed  and,  thus,  the  intensity  of  any 
beam  at  the  bottom  surface  of  the  crystal  can  be  determined,, 

The  theories  as  outlined  so  far  have  not  taken  into  account 

20 

absorption  or  relativistic  effects,  Hashimoto  et  alq  have  considered 
absorption  in  the  two-beam  case  by  assuming  a  complex  lattice  potential* 
Their  calculations  predict  that  the  only  effect  of  the  absorption  will 
be  a  decrease  in  intensity  but  no  change  in  extinction  distance „ 

There  have  been  no  published  calculations  on  the  effect  of  absorption 
on  extinction  distance  in  the  multi-beam  case.  However,  preliminary 
results  of  calculations  being  currently  carried  out  in  this  laboratory 
predict  that  when  only  two  strong  beams  are  present  any  corrections  to 
the  extinction  distance  would  be  second-order  in  magnitude.  Because 
many  electron  microscopes  have  accelerating  voltages  of  80  kV  or 
higher,  relativistic  effects  must  be  taken  into  consideration. 


■ 
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Fujiwara  developed  a  relativistic  dynamical  theory  of  electron  dif¬ 
fraction  using  the  Dirac  wave  equation*  His  results  show  that  the  non- 
relativistic  theory  developed  using  the  Schroedinger  wave  equation  can 
be  corrected  for  relativistic  effects  by  two  simple  substitutions. 

These  are  the  replacing  of  the  non-relativistic  wavelength  by  a  rela- 
tivistically  corrected  one  and  by  multiplying  the  Fourier  coefficient 
term  Ug  by  £  where 


e  =  (i 


Here  v  and  c  are  the  velocities  of  the  electrons  and  light  respectively. 


. 
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CHAPTER  III 

EXPERIMENTAL  TECHNIQUE 

Two  different  experiments  have  been  proposed  to  test  the 
predictions  of  the  two-beam  and  multi-beam  approximations  of  the  dyna¬ 
mical  theoryo  The  first  of  these  involved  the  measurement  of  changes 
in  the  spacing  of  extinction  contours  resulting  from  the  systematic 
variation  of  the  orientation  of  a  wedge-shaped  crystal.  From  such 
measurements  the  variation  of  $/£0  with  A9  can  be  determined.  The 
second  proposal  was  to  measure  the  actual  spacing  of  extinction  con¬ 
tours  in  the  images  of  wedge-shaped  crystals  of  known  wedge  angle. 

From  this  spacing,  the  magnification  of  the  image,  and  the  wedge  angle 
of  the  crystal,  values  of  6,Q  can  be  determined.  It  was  decided  to 
carry  out  these  experiments  using  the  {ill}  and  {220}  reflections  in 
silicon  at  electron  microscope  accelerating  voltages  of  100  and  152  kV. 

The  actual  procedure  in  making  these  measurements  can  be 
divided  into  four  parts »  These  are  calibration  measurements,  specimen 
preparation,  electron  microscope  examination,  and  analysis  of  observa¬ 
tions  . 

A.  Calibration  Measurements 


The  extinction  distance  measurements  were  carried  out  in  a 
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JEM  150  electron  microscope  using  a  tilting-rotating  stage  to  vary  the 
orientation  of  the  crystal  specimen  This  stage  allowed  the  specimen 
to  be  tilted  through  a  range  of  £  10°  about  an  axis  perpendicular  to 
the  incident  electron  beam  and  rotated  through  36O0  about  an  axis  per- 
pendicular  to  the  plane  of  the  specimen*  In  order  to  be  able  to  tilt 
the  specimen  through  a  known  angle s  a  high  precision  click- type  con¬ 
trol  that  allowed  the  angle  of  tilt  to  be  changed  in  increments  of 
0.025°  was  used. 

(i)  Tilt  Axis  Direction  Determination 

In  the  versus  A0  measurements,  the  specimen  was  oriented 

so  that  the  set  of  planes,  corresponding  to  a  reflection  of  interest, 
were  parallel  to  the  axis  of  tilt  of  the  stage*  In  this  orientation, 
a  given  change  in  the  angle  of  tilt  of  the  specimen  results  in  an  equal 
change  in  the  value  of  AO*  However,  in  order  to  orient  the  specimen  in 
this  manner,  the  direction  of  the  axis  of  tilt,  as  seen  at  the  viewing 
screen  of  the  microscope,  had  to  be  determined.  This  was  done  by  using 
the  fact  that,  as  the  specimen  is  tilted,  the  Kikuchi  pattern  must  move 
in  a  direction  perpendicular  to  the  tilting  axis  of  the  stage.  The 
direction  of  this  motion  was  found  by  taking  a  series  of  electron 
micrographs  of  the  diffraction  pattern  of  the  specimen  as  it  was  being 
tilted.  The  Kikuchi  patterns  in  these  micrographs  were  then  examined 
and  the  direction  of  motion  determined  from  the  direction  of  shift  of 


■ 
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the  patterns  in  the  micrographs •  This  method  of  determination  allowed 
the  direction  of  the  axis  of  tilt  as  seen  at  the  viewing  screen  to  be 
found  within  ±  3°®  It  should,  be  noted  that  this  direction  is  not  the 
same  as  the  actual  tilt  axis  direction  in  the  stage  due  to  rotation  of 
the  image  in  the  magnetic  lenses  of  the  microscope. 

(ii)  Tilt  Calibration 

The  precision  with  which  A.0  can  be  found  depends  on  the  ac¬ 
curacy  with  which  the  specimen  can  be  tilted.  The  specimen  tilting 
system  was  checked  for  both  long  range  errors,  which  would  result  in 
the  average  increment  of  tilt  being  different  than  the  designed  .025°» 
and  short  range  errors  which  would  result  in  random  variations  of  the 
increment.  The  long  range  errors  were  found  by  tilting  the  specimen 
through  a  known  angle  as  determined  from  the  initial  and  final  Kikuchi 
patterns  and  comparing  this  angle  with  the  angle  of  tilt  as  determined 
from  the  tilt  control. 

In  an  actual  check  the  specimen  was  oriented  so  that  a  row 
of  systematic  reflections  was  perpendicular  to  the  tilt  axis  direction. 
The  specimen  was  then  tilted  from  an  orientation  in  which  the  (hkl) 
reflection  was  in  the  Bragg  condition  to  one  in  which  the  (hkl)  re¬ 
flection  was  in  the  Bragg  condition  (see  Figure  5)«  The  exact  Bragg 
orientation  for  a  reflection  was  determined  from  the  fact  that,  under 
these  conditions,  a  bright  Kikuchi  line  must  pass  through  the  middle 
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of  the  spot.  Tilting  the  specimen  in  this  manner  results  in  the  ori¬ 
entation  of  the  specimen  being  changed  by  an  angle  equal  to  twice  the 
Bragg  angle  of  the  (hkl)  reflection.  The  angle  through  which  the 
specimen  was  tilted  as  given  by  the  tilt  control  was  then  compared  with 
this  angle  to  determine  the  accuracy  of  the  tilting  mechanism.  After 
carrying  out  this  check  over  different  parts  of  the  range  of  tilt  of 
the  stage,  it  was  found  that  the  angles  of  tilt  as  given  by  the  tilt 
control  were  within  one  per  cent  of  those  determined  from  Kikuchi 
patterns • 


The  short  range  random  variations  in  increment  size  were 
checked  by  using  a  method  given  by  Sheinin  •  In  general  they  were 
found  to  be  less  than  10  per  cent  of  the  increment  size  although  oc¬ 
casional  variations  ranging  up  to  50  per  cent  were  noted.  However, 
these  variations  usually  occurred  in  pairs  in  such  a  manner  that  the 
average  increment  found  from  the  two  of  them  would  still  be  approxi¬ 
mately  .025°*  Moreover,  such  variations  would  result  in  scatter  in 
the  results  of  a  size  which  is  already,  in  general,  less  than  the 
range  allowed  for  the  experimental  points  (see  Figures  10,  11,  12 
and  13)* 


The  accuracy  of  the  rotation  control  of  the  stage  was  not 


checked  as  precise  rotations  of  the  specimen  were  not  required 
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Figure  5 

A  drawing  of  a  diffraction  pattern  showing  the  initial  and  final  ori° 
entations  in  the  accuracy  check  of  the  stage#  The  arrow  shows  the 
direction  of  the  tilt  axis# 

(a)  Initial  orientation  showing  the  bright  (hkl)  Kikuchi  line  passing 
directly  through  the  (hkl)  spot# 

(b)  Final  orientation  showing  the  bright  (hkl)  Kikuchi  line  passing 
directly  through  the  (hkl)  spot 


#-.v  V- 


X  o 

(000) 


A(hkl) 

I  (b) 


FIGURE  5 


25 


(iii)  Magnification  Calibration 

In  measuring  £0  it  was  necessary  to  know  the  magnification 
of  the  images*  This  magnification  is  controlled  by  the  current  in  the 
intermediate  lens  of  the  microscope*  Thus,  in  order  to  know  the  mag¬ 
nification  of  the  image  it  was  necessary  to  accurately  calibrate  the 
magnification  of  the  microscope  against  the  intermediate  lens  current* 
This  was  carried  out  in  the  following  manner  using  a  carbon  replica  of 
a  diffraction  grating*  Electron  micrographs  of  the  replica  were  taken 
at  different  intermediate  lens  currents*  The  spacing  between  the  lines 
on  these  micrographs  was  then  measured  and.  compared  to  the  882  mp  per 
line  spacing  of  the  carbon  replica  in  order  to  determine  the  magnifica¬ 
tion*  A  series  of  micrographs  were  taken  at  each  intermediate  lens 
current  setting  in  order  to  minimize  the  effects  of  local  variations  in 
the  line  spacing  and  distortion  in  the  imaging  system*  This  calibration 
was  carried  out  at  both  100  and  152  kV  and  was  checked  once  during  the 
measurements  of  £0»  The  accuracy  of  this  calibration  was  limited  by  the 
tolerance  quoted  for  the  carbon  replica  and  the  precision  to  which  the 
intermediate  lens  current  could  be  determined*  These  two  combined  gave 
a  possible  error  of  +  3  per  cent  to  the  calibration  curve* 

B.  Specimen  Preparation 

In  order  to  make  precise  extinction  distance  measurements, 
smooth  regular  wedge-shaped  specimens  containing  no  strains  were 


- 
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required*  Silicon  was  chosen  as  the  material  from  which  to  prepare 
such  specimens  for  two  reasons*  Firstly,  a  technique  was  available  to 
grow  needle-shaped  crystals  or  whiskers  which  had  the  necessary  wedge 
shape  and,  secondly,  the  natural  brittleness  of  silicon  helped  to  pre¬ 
vent  the  formation  of  strains  in  these  specimens. 

The  technique  used  to  grow  the  whiskers  involved  the  dis¬ 
proportionation  of  Sil^.  The  actual  experimental  method  employed  was 
similar  to  that  used  by  other  researchers  ^3-25 #  jn  a  typical  run  the 
following  procedure  was  carried  out*  A  piece  of  clear  quartz  tubing 
(25  cm  in  length  and  about  140  ca?  in  volume)  was  sealed  at  one  end 
and  the  following  quantities  of  high  purity  chemicals  put  into  it? 

0,84  g  of  iodine,  1  g  of  silicon,  and  2  to  3  mg  of  nickel*  The  tube 
was  then  evacuated  to  a  pressure  of  about  5  x  10"^  Torr  and  sealed  off. 
In  order  to  prevent  the  loss  of  iodine  into  the  evacuating  system  at 
such  low  pressures,  the  end  of  the  tube  containing  the  chemicals  was 
immersed  in  liquid  nitrogen  during  the  evacuation  process.  The  sealed 
tube  was  then  put  in  a  horizontal  cylindrical  furnace  whose  operating 
temperature  varied  smoothly  from  1150°  C  at  one  end  to  800°  C  at  the 
other.  The  chemicals  in  the  tube  were  located  in  the  hotter  region  of 
the  furnace.  After  about  two  hours,  the  tube  was  withdrawn  and  allowed 
to  cool.  Inside  of  it  could  be  seen  whiskers  ranging  in  size  from 
about  a  centimeter  long  and  a  millimeter  in  diameter  down  to  very  fine 
filaments  less  than  a  micron  in  diameter.  The  size  cf  the  whiskers 


. 
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appeared  to  depend  on  the  temperature  of  the  region  in  which  they  grew 
with  the  larger  whiskers  being  located  in  the  hotter  regions.  After 
it  had  cooled,  the  tube  was  broken  open  and  the  whiskers  collected. 

26  27 

The  process  involved  in  the  growth  of  the  crystals  has  been  studied  *  ' 
and  found  to  be  a  thermo chemical  one  using  the  disproportionation  re¬ 
actions 

2SU2(g)  ^^Si(s)  +  Sil4(g) 

At  temperatures  in  the  range  of  1100°  C  the  reaction  moves  to  the  left 

while  at  temperatures  in  the  range  of  950°  C  the  reaction  moves  to  the 

right.  Thus,  in  the  furnace  this  reaction  results  in  a  net  transport 

of  silicon  from  the  hot  region  to  the  cooler  regions  of  the  tube. 

Under  the  conditions  described  above,  the  Si  comes  out  of  the  Sil2 
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compound  in  the  form  of  needle-shaped  crystals.  It  was  found  that 
the  addition  of  a  small  amount  of  impurity  such  as  gold  or  nickel 
greatly  increased  the  probability  of  formation  of  such  crystals. 

C.  Electron  Microscope  Examination 

The  determinations  of  extinction  distance  parameters  were 
made  from  observations  made  in  the  electron  microscope.  The  actual 
process  used  in  making  these  observations  consisted  of  choosing  and 
orienting  a  suitable  specimen  and  then  taking  a  series  of  electron 
micrographs  as  the  specimen  was  systematically  tilted.  However,  the 
choice  of  a  suitable  specimen,  how  it  oriented  and  how  the  micrographs 
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were  taken  depended  on  the  type  of  measurement  being  made, 

(i)  Measurements  of  S/£>0  Versus  A9 

A  specimen  suitable  for  making  this  type  of  measurement  was 
chosen  in  the  following  way,  First,  it  had  to  be  a  strain  free  region 
of  uniform  thin  wedge-shaped  single  crystal*  Thin  wedges  were  required 
so  that  the  extinction  contours  could  still  be  resolved  at  large  values 
of  A 0  where  £  is  quite  small®  Then,  a  diffraction  pattern  of  this 
region  had  to  contain  at  least  one  of  the  two  reflections  of  interest,. 
If  it  did,  the  specimen  was  then  oriented  so  that  the  set  of  planes 
corresponding  to  this  reflection  were  parallel  to  the  axis  of  tilt  of 
the  stage.  If  the  only  strong  low-order  reflections  then  present  in 
the  diffraction  pattern  were  in  the  row  of  systematic  reflections  cor¬ 
responding  to  this  set  of  planes,  the  specimen  was  used  for  making 
measurements  of  Zko9  This  last  feature  was  required  because  one  of 
the  purposes  of  this  experiment  was  to  test  the  predictions  of  the  two- 
beam  dynamical  theory.  The  closest  approximation  to  a  two-beam  situation 
achievable  in  a  real  crystal  is  that  in  which  only  a  single  row  of 
systematic  reflections  is  present  in  the  diffraction  pattern.  This  is 
difficult  to  achieve  in  practice  so  specimens  whose  diffraction  pat¬ 
terns  contained  high-order  spots  off  the  row  were  also  used.  Suitable 
regions  having  all  the  necessary  qualifications  were  usually  found  at 
the  intersection  of  two  whiskers  (see  Figure  6), 
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Figure  6 

An  electron  micrograph  showing  the  type  of  crystal  formation 
suitable  regions  for  <$/<S0  versus  A0  measurements  were  found. . 
filament  type  whiskers  are  also  visible,  x  500 
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The  variation  of  6/<$0  with  A0  was  found  from  micrographs  ob¬ 
tained  in  the  following  manner,,  After  a  suitable  region  of  specimen 
had  been  found  and  oriented  as  described  above,  it  was  tilted  until 
one  of  the  higher  order  systematic  reflections  was  in  the  Bragg  condi¬ 
tion,  as  determined  from  the  Kikuchi  pattern*  Using  this  orientation 
as  a  reference  point,  subsequent  orientations,  obtained  by  tilting  the 
specimen,  could  be  found  by  using  the  angle  of  tilt  determined  from 
the  tilt  control*  The  value  of  A0  for  a  parti cul.ar  orientation  was 
then  known  to  within  an  accuracy  of  £  0,025° »  the  sum  of  the  errors  in¬ 
volved  in  finding  the  reference  point  and  in  random  variations  of  the 
tilt  control.  From  the  reference  orientation,  the  specimen  was  then 
tilted  through  the  range  of  orientations  over  which  S/£0  versus  A0  was 
to  be  found.  In  this  range  it  was  tilted  in  steps  of  .025°  and  a  dark 
field  micrograph  using  the  reflection  of  interest  taken  at  each  step. 

In  order  to  be  able  to  find  Z/&0  from  these  micrographs,  the  range  of 
tilt  was  always  chosen  to  include  the  Bragg  orientation  of  the  reflec¬ 
tion  of  interest.  Also,  the  micrographs  were  taken  at  constant  magni¬ 
fication  and  with  the  condenser  lens  defocused  for  better  resolution. 
The  range  of  values  of  A0  over  which  these  micrographs  were  taken  was 
approximately  twice  the  Bragg  angle  of  the  reflection  of  interest 
centered  at  the  Bragg  condition  for  that  reflection.  For  each  series 
of  micrographs,  a  set  of  diffraction  patterns  through  the  same  range 
of  orientations  were  recorded  in  order  to  determine  the  presence  of  any 
low-order  strongly  diffracted  beams  other  than  those  in  the  row  under 
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consideration o 

(ii)  Measurements  of  £0 

These  measurements  require  that  the  wedge  angle  of  the  crys® 

tal  under  study  be  known.  Wagner  et  al?^  have  studied  the  morphology 

of  the  silicon  whiskers  produced  by  the  disproportionation  method  out® 

lined  above.  They  found  that  the  regular  whiskers  were  hexagonal  in 

cross  section  and  had  their  long  axis  in  a  [lllj  direction.  For 

whiskers  whose  diameters  were  less  than  a  micron,  the  lateral  faces 

were  fllO}  planes.  For  larger  whiskers,  however,  the  lateral  faces 

were  usually  {21l}  planes.  The  morphology  of  the  regular  whiskers 

grown  here  was  checked  by  comparing  micrographs  of  whiskers  with  their 

corresponding  electron  diffraction  patterns.  The  variation  in  thick® 

ness  across  the  whiskers  as  determined  from  the  micrographs  was  cor® 

related  with  the  diffraction  patterns  to  find  the  planes  forming  the 

lateral  faces.  The  results  of  these  measurements  for  submicron  size 
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whiskers  were  in  agreement  with  those  of  Wagner  et  al  « 

In  making  the  £0  measurements,  regular  whiskers  whose  ori® 
entation  and  dark  field  image  were  similar  to  that  shown  in  Figure  7 
were  used.  Because  of  the  large  wedge  angle  involved,  it  was  necessary 
to  use  only  whiskers  near  a  (110)  orientation  in  order  to  keep  rota¬ 
tions  of  the  whisker  about  its  long  axis  to  a  minimum.  At  least  three 
extinction  contours  were  required  in  the  image  of  a  wedge-shaped 
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Figure  7 

Configuration  of  whisker  and  image  in  S,Q  measurements, 

(a)  Perspective  drawing  of  a  regular  whisker  with  an  electron  beam 
striking  its  top  surface,  a  (110)  plane,  at  close  to  a  perpendicular 
angle • 

(b)  A  drawing  of  the  image  of  the  whisker  showing  the  presence  of 
extinction  contours  resulting  from  the  wedge-shaped  regions 
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region  for  accuracy  in  measuring  the  spacing. 

After  a  suitable  whisker  was  found.,  it  was  oriented  so  that 
one  of  the  (111)  or  (220)  reflections  present  in  a  (110)  diffraction 
pattern  was  near  its  Bragg  condition,  but  no  other  intense  low-order 
reflections  were  present,  A  dark  field  image  of  the  whisker  was  then 
obtained  and  its  magnification  found  from  the  intermediate  lens  cur¬ 
rent,  The  specimen  was  then  tilted  through  the  Bragg  orientation  in 
increments  of  equal  size  and  dark  field  micrographs  of  constant  magni¬ 
fication  taken  at  each  step,  A  series  of  diffraction  patterns  through 
the  same  range  of  orientations  were  also  recorded.  This  procedure  was 
then  repeated  for  the  other  reflection  of  interest, 

D.  Analysis  of  Observations 

In  order  to  find  both  the  variation  of  £/$0  with  A0  and  <g  , 
it  was  necessary  to  have,  at  least,  a  measure  of  the  spacing  between 
extinction  contours  in  the  dark  field  micrographs.  This  was  done 
through  use  of  a  microdensitometer  arrangement  (see  Figure  8),  In 
this  system  the  micrographs  were  carried  at  constant  speed  through  a 
light  beam  and  the  variations  in  transmitted  intensity  converted  into 
peaks  on  a  chart  recorder  (see  Figure  9)*  The  spacing  was  then  found 
by  measuring  the  distance  between  the  peaks.  However,  before  a  set  of 
traces  were  taken,  the  micrographs  were  examined  for  areas  where  the 
contours  were  evenly  spaced  and  where  there  was  no  evidence  of  strains. 
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Figure  8 


Microdensitometer  system  for  measuring  the  distance  between  extinction 


contours 
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The  position  of  such  an  area  was  marked  on  all  the  micrographs  in  the 
set#  These  micrographs  were  then  placed  in  the  microdensitometer  so 
that  the  contours  were  perpendicular  to  the  direction  of  travel  and 
traces  taken  of  the  marked  areas# 

(i)  Measurements  of  f/4  Versus  A0 

The  actual  variation  of  £/£0  with  A0  was  determined  from  the 
average  distance  between  the  peaks  in  the  intensity  profiles#  This 
distance  was  found  from  measurements  carried  out  over  a  limited  region 
of  the  traces  in  order  to  minimize  the  effects  of  uneven  thickness  var¬ 
iations  (see  Figure  9)»  The  average  spacing  of  the  peaks  at  the  Bragg 
condition  was  found  by  an  interpolation  process  involving  the  averages 
measured  in  the  vicinity  of  this  orientation#  This  average  at  the 
Bragg  condition  was  then  divided  into  the  measured  average  from  each 
trace  to  give  its  corresponding  value  of  These  values  were  then 

plotted  up  and  the  position  of  the  exact  Bragg  orientation  on  this 
curve  was  found  by  combining  the  information  given  by  the  initial  ori¬ 
entation  or  reference  point  with  the  symmetry  of  the  curve# 

(ii)  Measurements  of  £0 

In  the  measurements,  it  was  not  necessary  to  know  the 

actual  distance  between  the  contours  on  the  micrographs#  However,  in 
measurements  of  £0  this  spacing  had  to  be  found#  This  was  done  from 
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Figure  9 

Microdensitometer  traces  showing  the  spacing  of  the  extinction  contours 

(a)  near  the  Bragg  condition® 

(b)  at  a  A0  of  0.275°. 

Measurements  of  the  average  peak  to  peak  distance  would  be  carried  out 
in  the  region  A  to  B  on  all  traces  corresponding  to  this  set  of  micro¬ 
graphs 
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the  densitometer  traces  by  taking  into  account  the  rates  of  travel  of 
both  the  micrograph  through  the  light  beam  and  the  chart  paper  in  the 
recorder.  for  each  trace  was  then  calculated  from  this  distance  by 
taking  into  consideration  the  wedge  angle  and  rnagnifi cation *  The 
value  of  d0  corresponding  to  the  set  of  traces  was  then  found  from  the 
values  of  by  an  interpolation  process. 
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CHAPTER  IV 
RESULTS 

The  extinction  distance  measurements  were  carried  out  a  num¬ 
ber  of  times  for  each  reflection  at  each  accelerating  voltage.  The 
predictions  of  the  two-beam  and  multi-beam  dynamical  theories  of  elec¬ 
tron  diffraction  were  then  compared  with  the  averages  of  these  measure¬ 
ments# 

A#  Results  of  £/<50  Versus  A0  Measurements 

The  results  of  these  measurements  are  presented  in  Figures 
10 »  11,  12  and  13  along  with  the  corresponding  predictions  of  the  two- 
beam  and  multi-beam  theories.  The  experimental  results  presented  in 
these  Figures  are  the  averages  of  3»  4  and  5  sets  of  measurements 

respectively.  Only  those  results  were  used  in  which  the  variation  of 
£/£0  with  A0  was  relatively  smooth,  showing  no  recognizable  effects 
resulting  from  non- systematic  reflections.  These  effects  could  be 
detected  by  the  sharp  changes  which  occurred  in  the  £/<S0  versus  A0 
curves  (see  Figure  14),  These  changes  were  found  to  be  associated 
with  a  non-systematic  reflection  near  its  Bragg  condition*  The  exper¬ 
imental  averages  are  given  in  the  Figures  10  to  13  with  a  range  of 
error  of  £  2-J  per  cent#  The  experimental  results  all  fell  within  this 
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Figure  10 


The  variation  of  £/<£0  with  A0  for  the  (220)  reflection  in  silicon  at 


100  k V 
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Figure  11 


The  variation  of  €,!€>Q  with  A0  for  the  (220)  reflection  in  silicon  at 
152  kV 
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Figure  12 

The  variation  of  $/<$0  with  A0  for  the  (111)  reflection  in  silicon  at 
100  kV*  The  sharp  drops  in  <$/<50  take  place  at  approximately  +  0,260° 
and  -  0,300° 
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Figure  13 


The  variation  of  £/£0  with  A.0  for  the  (111)  reflection  in  silicon  at 
152  kV*  The  sharp  drops  in  £/£0  take  place  at  approximately  ±  0#210C 


' 


o 


°5/5 


o 

o 

o 

o 

o 

o 

o 

o 

00 

o 

wo 

o' 

o 

o 

o 

o 

A 9  (DEGREES) 
FIGURE  13 


43 


Figure  14 

A  set  of  (111)  experimental  results  showing  the  effects  of  a  non- 
systematic  reflection  in  the  region  of  A0  from  0*050°  to  0.150° 
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range  except  where  shown  by  the  extended  error  bars.  This  range  of 
error  of  £  Zj  per  cent  arises  in  the  measurement  of  the  distance  between 
peaks  in  the  densitometer  traces  due  to  the  difficulty  in  defining  the 
exact  position  of  the  peaks. 

The  predictions  of  the  two -be am  and  multi-beam  theories, 
given  in  Figures  10  to  13,  were  calculated  using  the  expressions  devel¬ 
oped  earlier.  The  scattering  factors  used  in  the  calculations  were 

28 

found  in  tables  compiled  by  Ibers  and  Vainshtein  •  In  making  the 
multi-beam  calculations,  only  systematic  reflections  were  considered. 

The  number  of  these  reflections  used  in  a  calculation  was  determined 
by  a  trial  process  in  which  calculations  were  made  with  an  increasing 
number  of  beams,  until  the  addition  of  more  beams  changed  the  extinc¬ 
tion  distance  by  less  than  1  per  cent.  In  the  (111)  multi-beam  calcu¬ 
lations  this  condition  required  the  use  of  systematic  reflections  from 
(333)  to  (666)  while  in  the  (220)  calculations  systematic  reflections 
from  (440)  to  (880)  were  required. 

The  extinction  distances  in  the  multi-beam  calculations  were 
found  from  intensity  profiles  (see  Figure  15) •  These  are  plots  of  the 
variation  of  the  intensity  of  the  beam  under  consideration  with  depth 
in  the  crystal,  as  obtained  from  an  autoplotter  system  used  in  conjunc¬ 
tion  with  a  computer.  The  extinction  distance  was  found  by  measuring 
the  average  distance  between  the  peaks.  Attenuation  was  not  taken  into 
consideration  in  these  calculations  as  preliminary  calculations  showed 
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Figure  15 


An  intensity  profile  of  the  variation  with  depth  of  the  (111)  dif¬ 
fracted  beam  intensity  at  152  kV  as  found  from  a  10-beam  approximation 
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that  it  would  result  in  less  than  a  1  per  cent  correction  to  the  cal¬ 
culated  extinction  distances* 

The  actual  effect  of  the  systematic  reflections  on  the  extinc¬ 
tion  distance  is  different  for  the  (220)  and  (111)  reflections*  As 
can  be  seen  in  Figures  10  and  11,  the  predictions  of  the  two-beam  and 
multi-beam  theories  in  the  case  of  the  (220)  reflections  are  very  simi¬ 
lar*  Moreover,  the  results  of  the  experimental  measurements  are  in 
good  agreement  with  these  predictions*  However,  in  the  case  of  the 
(111)  reflection,  the  two-beam  and  multi-beam  theories  show  a  marked 
difference  in  the  dependence  of  <5/£0  on  AO*  The  two-beam  theory  pre¬ 
dicts  sharp  drops  in  the  extinction  distance  at  angular  deviations  of 
approximately  jh  0*280  and  £  0*225  degrees  at  100  and  152  kV  respec¬ 
tively*  In  the  regions  where  the  drops  take  place,  the  variation  of 
diffracted  beam  intensity  with  depth  becomes  complex  and  the  peaks  are 
not  clearly  defined  ( see  Figure  16) »  Thus  a  determination  of  extinc¬ 
tion  distance  is  ambiguous*  The  results  of  the  (111)  measurements 
show  that  the  experimental  extinction  distances  have  drops,  occurring 
in  the  same  general  range  of  A0*s,  which  are  similar  to  the  ones  pre¬ 
dicted  by  the  multi -beam  theory*  Moreover,  in  the  regions  where  the 
drops  occur,  the  intensity  variation  becomes  complex  (see  Figure  17)* 
However,  the  experimental  extinction  distances  at  both  100  and  152  kV 
decrease  more  rapidly  than  either  the  two-beam  or  multi-beam  theories 
predict*  The  results  at  100  kV  show  some  asymmetry  with  respect  to 
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Figure  16 

Complex  depth  periodicity  of  the  (111)  reflection  at  152  kV  as  pre¬ 
dicted  by  the  10-beam  approximation  for  A0  =  0*250 


Figure  17 

A  densitometer  trace  of  the  extinction  contours  corresponding  to  the 
(111)  reflection  at  152  kV  showing  the  complex  periodicity  at 
AO  a  0*200° 
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this  decrease  in  that  the  extinction  distance  appears  to  decrease  more 
rapidly  for  positive  values  of  A0  than  for  negative  ones. 

B.  Results  of  the  £0  Measurements 

Z0  measurements  were  attempted  for  the  (111)  and  (220)  re- 
flections  at  both  100  kV  and  152  kV«  However,  consistent  results  were 
only  obtained  for  the  (111)  reflection  at  152  kV«  Here,  the  measure¬ 
ments  gave  a  value  of  840  £  40  A  for  <?0  as  compared  to  the  two-beam 
theory* s  value  of  720  A  and  the  10-beam  approximation's  value  of  736  A. 
Results  were  also  obtained  for  the  (220)  reflections  at  152  kV  but  due 
to  multi-beam  effects  the  scatter  was  quite  large •  The  results  did 
indicate,  however,  that  this  extinction  distance  was  in  fairly  good 
agreement  with  the  two-beam  value  of  880  A  and  the  7“beam  value  of 
840  A •  No  results  were  obtained  at  100  kV  for  either  reflection  due 
to  multi-beam  effects.,  These  effects  made  the  periodicities  so  com¬ 
plex  and  of  such  low  contrast  that  the  defining  of  an  extinction  dis¬ 
tance  was  questionable 0 
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CHAPTER  V 
DISCUSSION 


The  results  of  the  extinction  distance  measurements  indicate 
both  agreement  and  disagreement  with  theory,  depending  on  which  re» 
flection  is  being  considered.  The  (111)  extinction  distance  varies 
with  A0  in  a  manner  different  from  the  predictions  of  both  the  two** 
beam  and  multi-beam  approximations,  €>  measurements  for  this  reflec¬ 
tion  also  disagree  with  both  approximations.  The  (220)  extinction 
distance,  however,  shows  good  agreement  with  the  two  approximations 
for  both  versus  A0  and  £  , 

A,  (111)  Extinction  Distance  Measurements 

As  seen  in  Figures  12  and  13,  the  (ill)  experimental  extinc¬ 
tion  distance  showed  a  sharp  decrease  similar  to  that  predicted  by 
the  10-beam  approximation.  However,  the  size  of  the  drop  was  greater 
than  that  predicted  by  this  theory.  Also,  the  measured  value  of  £Q 
for  this  reflection  was  20  per  cent  larger  than  the  theory  predicted, 

A  possible  explanation  of  these  results  lies  in  the  values 

P  ft 

of  the  scattering  factors^  used  in  making  the  theoretical  calcula- 
29 

tions.  Ibers  has  given  a  discussion  of  the  calculation  of  the  scat¬ 
tering  factors  and  the  sources  of  error  involved.  For  silicon,  the 
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factors  were  calculated  using  Hartree-Fock  self-consistent  wave  func¬ 
tions,  In  using  this  type  of  function,  he  pointed  out  that  there  were 
two  possible  sources  of  error.  The  first  of  these  were  errors,  either 
of  omission  or  of  calculation,  in  the  contributions  to  the  atomic 
fields  of  the  outer  electrons.  The  second  were  errors  arising  from 
the  use  of  Hartree-Fock  wave  functions,  calculated  for  spherically 
neutral  atoms,  for  an  atom  in  a  crystal  lattice.  In  a  later  paper^ 
he  gave  the  following  limits  to  the  accuracies  of  the  scattering  fac¬ 
tors  due  to  these  errors.  For  values  of  sin  0/A  greater  than  0,3  A“^ 
the  factors  are  reliable  to  10  per  cent  but  for  sin  0/A  less  than 
0,1  the  error  may  be  as  large  as  50  to  100  per  cent.  Here,  0  is 
one  half  of  the  scattering  angle  and  A  is  the  wavelength  of  the  inci¬ 
dent  electron.  In  the  case  of  the  (111)  reflection  in  silicon  sin  0/A 
is  equal  to  0,159  A~"*".  Thus,  the  scattering  factor  for  this  re¬ 
flection  can  only  be  considered  accurate  to  approximately  30  per  cent. 

An  experimental  value  for  the  (111)  scattering  factor  was 
calculated  by  varying  the  value  of  the  (111)  scattering  factor  used 
in  the  10-beam  approximation  until  a  value  of  840  A  was  found  for  CQ, 
the  (111)  extinction  distance  at  152  kV,  The  value  of  scattering  fac¬ 
tor  found  in  this  way  was  2,83  A,  This  experimental  value  is  within 
14  per  cent  of  the  value  of  3„26  A  given  by  Ibers  and  Vainshtein, 

Using  this  new  value  of  scattering  factor  for  the  (111)  reflection, 
the  variation  of  £  with  A0  was  recalculated  in  the  case  of  the  10-beam 
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approximation  and  compared  with  the  experimental  results  (see  Figures 
18  and  19)®  The  good  agreement  between  the  experimental  results  and 
the  new  empirical  curve  is  especially  noticabl©  in  Figure  19  which 
shows  the  results  at  152  kV®  As  was  seen  in  Figure  13,  the  experi¬ 
mental  results  dropped  below  both  the  two-beam  and  10-beam  curves  even 
at  angular  deviations  smaller  than  that  at  which  the  sharp  drops  took 
place®  However,  as  seen  in  Figure  19 »  the  new  empirical  curve  passes 
through  all  the  experimental  points  between  -0«175  and  +  0 « 175 ®  More¬ 
over,  the  size  of  the  sharp  drop  as  predicted  by  the  empirical  curve 
gives  better  agreement  than  the  theoretical  curves  of  Figure  13®  Thus, 
it  can  be  seen  that  the  C/£0  versus  A.0  measurements  confirm  an  experi¬ 
mental  value  of  the  scattering  factor  of  2. ,83  A® 

At  100  kV  the  fit  between  the  experimental  results  and  em¬ 
pirical  curve  is  not  as  good  as  at  152  kV  (see  Figures  18  and  19)® 

Here,  the  asymmetry  of  the  experimental  results  at  100  kV  is  again 
demonstrated  as  the  results  at  negative  A0,  before  the  drop  takes 
place,  fall  above  the  empirical  curve,  while  for  the  corresponding 
values  of  positive  A0,  the  results  fall  on  the  curve®  A  possible  ex¬ 
planation  for  this  asymmetry  lies  in  the  effect  of  non-systematic 
reflections  arising  from  the  orientations  of  the  specimens®  In  the 
specimens  chosen  for  this  set  of  £/$0  measurements,  it  was  observed 
that  low-order  non-systematic  reflections  approached  and  passed 
through  their  Bragg  conditions  as  the  specimens  were  tilted  in  the 
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Figure  18 

A  comparison  of  the  variation  of  £/<50  with  A0  as  measured  for  the  (111) 
reflection  at  100  kV  with  the  corresponding  predictions  of  a  10-beam 
approximation  calculated,  using  the  experimental  value  of  2.83  A  for  the 
(111)  scattering  factor 
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Figure  19 


A  comparison  of  the  variation  of  £/£>0  with  A0  as  measured  for  the  (111) 
reflection  at  152kV  with  the  corresponding  predictions  of  a  10-beam 


approximation  calculated  using  the  experimental  value  of  2,83  A  for 
the  (111)  scattering  factor 
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direction  of  negative  A0.  A  similar  effect  was  not  noticed  at  152  kV, 
probably  due  to  the  different  curvature  of  the  Ewald  sphere.  Also,  it 
was  not  observed  in  the  case  of  the  (200)  reflection  where  different 
orientations  were  involved.  Thus,  the  results  at  100  kV  are  not  incon¬ 
sistent  with  the  experiment  scattering  factor  of  2,83  A. 


The  use  of  the  experimental  scattering  factor  of  2.83  A  is 
also  supported,  by  independent  measurement s^“33  carried  out  to  find  the 
corresponding  X-ray  scattering  factors.  The  electron  scattering  factor, 
fe,  is  related,  to  the  X-ray  scattering  factor,  fx,  through  the  rela- 
tion28 


0.023934 - £— 

(sin  9)2/At 


where  Z  is  the  atomic  number  of  the  element  being  considered  and  0  and 

34 

X  are  as  defined  above.  The  most  reliable  of  the  results  are  those 

33 

of  Hattori  gt  &1:  who  found  the  fx  factors  for  different  reflections 

in  silicon  from  the  spacing  of  Pendellosung  fringes,  a  phenomenon  in 

X-ray  diffraction  similar  to  that  of  extinction  contours  in  electron 

34 

diffraction.  Their  results,  as  corrected  by  Dawson-'  ,  gave  a  value  of 
10.90  £  .04  for  the  fx  of  the  (111)  reflection  as  compared  to  Ibers' 
and  Vainshtein's  value  of  10,54.  Substituting  this  value  into  the 
expression  for  fQ,  a  value  of  2,92  A  is  obtained  for  the  electron  scat¬ 
tering  factor.  This  agrees  quite  favourably  with  the  empirical  value 
of  2,83  A. 
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B»  (220)  Extinction  Distance  Measurements 

The  predictions  of  the  two-beam  and  7-beam  approximations 
are  so  similar  for  this  reflection  that  it  is  impossible,  from  the 
experimental  results,  to  determine  which  approximation  best  describes 
the  actual  physical  situation*  The  experimental  results  agree  well 
with  both  approximations  at  both  100  and  152  kV»  The  lack  of  consist¬ 
ent  results  for  the  measurements  of  £0  for  this  reflection  prevents 
correlation  of  these  with  the  S  versus  A0  measurements*  However, 
Frankl  ,  in  his  measurements  of  the  (220)  extinction  distance  in  sil¬ 
icon  found  a  value  of  ?20  A  at  100  kV  for  £>q9  in  excellent  agreement 
with  that  predicted  by  the  7-beam  approximation  of  723  A*  Thomas  and 
Levine-^  also  made  measurements  of  the  (220)  extinction  distance  using 
stacking  fault  fringes  rather  than  thickness  contours*  They  obtained 
a  value  of  660  A  for  £  at  100  kV®  They  did  not  explain  the  deviation 
from  the  theoretical  value  although  they  did  say  that  it  was  -under 
further  investigation*  Booker  and  Hazzledine-^  have  investigated  the 
profiles  of  planar  defects  in  crystals  such  as  stacking  faults®  Due 
to  the  different  variables  such  as  the  defect's  phase  angle,  the  angu¬ 
lar  deviation  of  the  reflection  being  considered  from  its  Bragg  con¬ 
dition,  and  absorption  of  the  electron  beams  in  the  crystal,  they 
concluded  that  "the  greatest  care  must  be  taken  in  interpreting  elec¬ 
tron  micrographs  of  such  faults  before  definite  conclusions  can  be 
reached®"  Thus,  the  more  reliable  results  are  those  of  Frankl  which 


56 


are  also  consistent  with  the  £  versus  40  measurements  made  here. 

The  X-ray  scattering  factor  has  also  been  measured  for  this 
31-33 

reflection^  ^  and  found  to  vary  slightly  from  that  given  by  Ibers 
28 

and  Vainshtein  •  However,  this  variation  results  in  only  a  3'i'  per 
cent  change  in  the  electron  scattering  factor.  This  change  in  scatter¬ 
ing  factor  is  too  small  to  be  detected  here  because  it  lies  within  the 
range  of  error  of  the  present  experiments. 

C.  Multi -Be am  Effects 

The  effects  of  beams  other  than  the  two  principal  beams  on 
the  extinction  distance  can  be  considered  for  the  two  cases  of  system¬ 
atic  and  non-systematic  reflections. 

The  effects  of  systematic  reflections  were  taken  into  con¬ 
sideration  in  the  multi -beam  calculations.  In  the  case  of  the  (220) 
reflection,  the  7-beam  calculations  predicted  values  of  £0  of  72 3  and 
840  Angstroms  at  100  and  152  kilovolts  respectively  as  compared  to  the 
two-beam  values  of  760  and  890  Angstroms.  However,  the  multi-beam 
predictions  of  the  variation  of  S,  with  40  differed  from  that  predicted 
by  the  two-beam  theory  by  an  amount  too  small  to  detect  by  the  meas¬ 
urements  carried  out  here.  Thus,  for  the  (220)  reflection  the  effect 
of  the  systematic  reflections  is  to  decrease  but  to  have  a  rela¬ 
tively  small  effect  on  the  variation  of  £  with  40. 
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In  the  case  of  the  (111)  reflection  the  effects  of  systematic 

reflections  were  important  for  both  ZQ  and  £  versus  A 0.  For  £  the 

10 -be am  approximation  predicted  values  of  63O  and  736  Angstroms  at  100 

and  132  kilovolts  respectively  as  compared  to  the  two-beam  values  of 

615  and  720  Angstroms.  Thus,  for  the  (111)  reflection  the  effect  of 

systematic  reflections  is  a  predicted  increase  in  <|0  as  opposed  to  the 

usual  decrease-^.  In  the  £  versus  A0  measurements  the  two-beam  theory 

predicts  a  smooth  continuous  variation  of  £  with  A0,  with  no  sharp 

drops.  However,  3"beam  calculations  which  included  either  the  (Ill) 

or  (222)  reflection  were  carried  out  and  in  these  cases  sharp  drops 

were  predicted.  Moreover,  these  drops  took  place  in  the  same  range 

of  values  of  A0  as  they  were  observed  experimentally.  As  was  shown 

earlier  in  Figures  12  and  13,  when  10  beams  were  considered  these 

drops  were  also  predicted.  However,  the  3”6eam  approximation  gave 

qualitative  agreement  with  the  experiment  results  and  this  agreement 

was  only  improved  slightly  by  using  a  10-beam  approximation.  Thus, 

this  effect  is  probably  associated  with  the  " forbidden"  (222)  reflec- 

39 

tion  in  diamond  structure^  as  this  reflection  first  appears  in  the 
off-diagonal  elements  of  the  3“team  approximations.  This  reflection 
is  called  "forbidden,"  even  though  it  is  observed  in  electron  dif¬ 
fraction  patterns,  because  its  structure  factor  is  zero.  Heidenreich' 
explained  the  presence  of  a  (222)  reflection  by  double  diffraction 
involving  a  pair  of  reflections  such  as  (113)  and  (111).  However, 

(222)  reflections  were  observed  to  be  present  even  when  pairs  of  re- 
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flections  such  as  these,  which  could  combine  to  give  a  (222)  reflection, 

40  41 

were  not*  Hoerni  and.  Fujimoto  explained  the  presence  of  the  (222) 
reflection  by  taking  into  account  interactions  between  the  systematic 
reflections  including  the  (111)  reflection.*  These  interactions  were 
calculated  to  contribute  the  main  part  of  the  (222)  intensity  even  when 
diffracting  pairs  such  as  (113)  and  (111)  were  present,  Thus,  the 
presence  of  the  ''forbidden'1  (222)  reflection  is  predicted  to  affect  the 
intensity  of  the  (111)  reflection,  possibly  through  sharp  changes  in 
extinction  distance  as  were  observed® 

The  effects  of  non-systematic  reflections  are  much  more  dif¬ 
ficult  to  take  into  consideration  as  there  is  no  accurate  method  of 
determining  how  far  these  reflections  are  from  their  exact  Bragg  con¬ 
ditions  for  any  random  orientation®  As  was  seen  in  Figure  14,  one  of 
the  possible  effects  of  these  reflections  is  to  change  the  variation 
of  S,  with  A.9  by  introducing  small  sharp  changes  in  extinction  distance® 
Other  effects  were  observed  in  making  the  measurements.  Here,  the 
non-systematic  reflections  were  found  to  both  reduce  the  contrast  in 
the  image,  thus,  making  it  difficult  to  differentiate  between  the 
bright  and  dark  contours,  and  to  introduce  scatter  into  the  measured 

results  particularly  in  the  case  of  the  (220)  reflections  at  152  kV, 

42 

These  observations  are  supported  by  work  done  by  Humphreys  et  al. 
on  stacking  fault  fringes  in  gold  and  silicon.  They  carried  out  cal¬ 
culations  and  measurements  of  the  effects  of  both  systematic  and  non- 
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systematic  reflections  and  obtained  results  similar  to  those  given 
here# 

D.  Conclusions 

The  original  purpose  of  the  work  undertaken  here  was  to  find 
the  variation  of  extinction  distance  with  angular  deviation  from  the 
Bragg  condition  for  the  (111)  and  (220)  reflections  in  silicon  and  to 
compare  the  results  with  the  predictions  of  the  two-beam  and  multi¬ 
beam  approximations  of  the  dynamical  theory  of  electron  diffraction® 

In  the  case  of  the  (111)  reflection,  a  10-beam  approximation 
involving  only  systematic  reflections  and  calculated  using  2® 83  A  for 
the  (111)  electron  scattering  factor,  a  value  obtained  from  %Q 
measurements ,  was  found  to  be  in  good  agreement  with  the  experimental 
results.  The  use  of  this  experimental  scattering  factor  rather  than 
the  theoretical  value  of  3»26  A  is  supported  by  recent  measurements 
made  of  the  X-ray  scattering  factor  for  the  same  reflections#  The 
measured  X-ray  factor  corresponds  to  a  value  for  the  electron  scatter¬ 
ing  factor  of  2®92  A  which  is  in  close  agreement  with  the  experimental 
value  found  here.  The  10-beam  approximation  also  predicts  the  observed 
sharp  drops  in  extinction  distance,  a  phenomenon  probably  related  to 
the  ’’forbidden"  (222)  reflection®  In  the  case  of  the  two-beam  approxi¬ 
mation  there  is  poor  agreement  with  the  experimental  results® 
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In  the  ease  of  the  (220)  reflection  the  predictions  of  the 
two-beam  and  7-beam  approximation  differ  only  slightly  and  are  both  in 
good  agreement  with  the  experimental  results*  These  £  versus  A0  meas¬ 
urements  were  also  in  agreement  with  the  £0  measurements  of  Frank! 0 

Other  workers  have  found  similar  results  to  those  obtained  here*  Uyeda 
12 

and  Nonoyama  found  agreement  between  the  two-beam  theory  and  experi¬ 
mental  results  for  the  (200)  reflection  in  MgO.  Within  the  range  of 

43 

angular  deviations  considered  here,  Sheinin  also  found  good  agreement 
between  experimental  results  and  the  predictions  of  the  two-beam  and 
12-beam  approximations  for  the  (110)  extinction  distance  in  molybdenum* 
Thus,  the  variation  of  €  with  A  6  for  the  (111)  reflection  is  unique 
among  published  measurements  of  this  type  in  that  it  differs  markably 
from  the  two-beam  predictions*  As  pointed  out  earlier,  this  variation 
is  probably  due  to  the  "’forbidden"  (222)  reflection* 

Eo  Suggestions  For  Further  Work 

More  measurements  of  the  variation  of  extinction  distance 
with  deviation  from  the  Bragg  condition  for  different  reflections  in 
silicon  and  similar  measurements  in  other  materials  might  be  carried 
out  to  check  the  effect  of  forbidden  reflections*  However,  it  should 
be  noted  that  extinction  distance  is  basically  a  parameter  of  the  two- 
beam  theory  arising  out  of  the  predicted  sine  squared  variation  for 
the  intensity  of  a  beam  with  depth  in  a  crystal.  The  presence  of  other 
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reflections  affects  both  the  shape  of  the  intensity  peaks  and  the  spac¬ 
ing  between  them*  The  severity  of  these  effects  depends  on  both  the 
proximity  of  the  reflections  to  their  Bragg  conditions  and  the  energy 
of  the  incident  electrons*  The  latter  condition  arises  from  the  fact 
that  systematic  reflections  would  be  expected,  to  have  a  greater  effect 
as  the  energy  of  the  electron  beam  increases'^,  Thus,  the  use  of  a 
parameter  such  as  extinction  distance  to  describe  a  multi-beam  case  is 
open  to  question*  Perhaps,  a  more  accurate  method  of  checking  the 
various  theories  is  to  compare  the  shape  of  the  theoretical  depth 
intensity  profiles  of  a  beam  with  the  actual  observed  shape  in  the 
densitometer  traces.  This  approach  would  be  especially  useful  in  the 
regions  of  complex  periodicity*  However,  in  making  theoretical  pre¬ 
dictions  of  the  shape  of  contours  it  is  necessary  to  include  the 
effects  of  absorption*  Such  effects  are  not  yet  completely  understood. 
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